Astronomy & Astrophysics manuscript no. paper 


February 2, 2008 


(DOI: will be inserted by hand later) 





On the progenitor of V838 Monocerotis 

R. Tylenda^''^, N. Soker^, and R. Szczerba^ 



in 
o 
o 

(N 



Department for Astrophysics, N.Copernicus Astronomical Center, Rabiariska 8, 87-100 Toruri, Poland 
e-mail: tylenda@ncac.torun.pl, szczerba@ncac.torun.pl 
^ Department of Physics, Technion-Israel Institute of Technology, 32000 Haifa, Israel 
e-mail: sokerSphysics . technion. ac . il 

Centre for Astronomy, N. Copernicus University, 87-100 Torun, Poland 
Received 



m 
> 
m 
oo 

(N 
O 



Abstract. We summarize and analyze the available observational data on the progenitor and the enviroment of 
V838 Mon. From the available photometric data for the progenitor of V838 Mon we exclude the possibility that 
the object before eruption was an evolved red giant star (AGB or RGB star). We find that most likely it was a 
main sequence or pre-main sequence star of ~ 5 — lOM© . From the light echo structure and evolution we conclude 
that the reflecting dust is of interstellar nature rather than blown by V838 Mon in the past. We discuss the 
IRAS and CO data for interstellar medium observed near the position of V838 Mon. Several interstellar molecular 
regions have radial velocities similar to that of V838 Mon, so dust seen in the light echo might be related to one 
of them. 
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1. Introduction 

V838 Mon is a star caught in eruption at the beginning 
of January 2002 (Brown 1200211 ■ The eruption, as observed 
in optical wavelengths, lasted about three months, and 
was composed of two or three major peaks. After develop- 
ing an A-F supergiant spectrum at the optical maximum 
at the beginning of February, the object showed a gen- 
eral tendency to evolve to lower effective temperatures. In 
April 2002 it almost disappeared from the optical but re- 
mained very bright in infrared becoming one of the coolest 
M-type supergiants yet observed. Detailed descriptions of 
the spectral and photometric evolution of V838 Mon can 
be found in a number of papers including Munari et al. 
H2002b|l . Kimeswenger et al. (2002), Kolev et al. PUH^ . 
Osiwala et al. II2003|I . Wisniewski et al. H2003|l . Crause 
et al. ((2003) and Kipper et al. (2UU^. 

The nature of the V838 Mon eruption is enigmatic. 
As discussed by Soker & Tylenda (|2003|l . thermonuclear 
models (classical nova, He-shell flash) seem to be unable to 
explain this type of eruption. Therefore other mechanisms 
such as a stellar merger model (Soker & Tylenda 2003 ) or a 
giant swallowing planets scenario (Retter & Marom .2003 ) 
have been proposed. 
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The global fading of V838 Mon in optical after out- 
burst has enabled us to discover a faint hot continuum in 
short wavelengths (Desidera & Munari 120021 Wagner & 
Starrfield mUT) later classified as coming from a normal 
B3 V star (Munari et al. I2002ajl . This strongly suggests 
that V838 Mon is a binary system which can be an im- 
portant fact for identifying the outburst mechanism. 

V838 Mon has received significant publicity due to its 
light echo, which was discovered shortly after the main 
eruption in February 2002 (Hcnden et al. I2002|) . and was 
seen in images by the HST (Bond et al. lM75|l . The hght 
echo was used, e.g. in Bond et al. H2003I OHinjl . to claim 
that the echoing matter was ejected by V838 Mon in pre- 
vious eruptions. This conclusion was disputed by Tylenda 
()2004|) , who examined the evolution of the light echo and 
concluded that the dust illuminated by the light echo was 
of interstellar origin rather than produced by mass loss 
from V838 Mon in the past. 

van Loon et al. (2004) argue that there are mul- 
tiple shells around V838 Mon, which were ejected by 
V838 Mon in previous eruptions. Hence they reason that 
prior to eruption V838 Mon was an asymptotic giant 
branch (AGB) star. These authors have also analyzed the 
light echo with more recent observations than in Tylenda 
f2004'), and argue that the echoing dust was ejected by 
V838 Mon in past eruptions. 
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Table 1. Photometry of the V838 Mon progenitor 



Band 


Xef f (/ini) 


Magnitude 


Reference 


B 


0.44 


15.87±0.10 


Kimesweneer et al. 120021 


B 


0.44 


15.81±0.06 


Goranskii et al. 120041 


Rc 


0.65 


14.84±0.06 


Goranskii et al. 120041 


R 


0.71 


14.56±0.10 


Kimeswenffer et al. 120021 


Ic 


0.80 


14.27±0.03 


Goranskij et al. PIMI 




0.82 


14.51±0.03 


DENIS 


J 


1.25 


13.86±0.11 


DENIS 


J 


1.25 


13.87±0.05 


2MASS 


H 


1.65 


13.51±0.06 


2MASS 


Ks 


2.15 


13.14±0.16 


DENIS 


Ks 


2.17 


13.33±0.06 


2MASS 



the bands for which two measurements have been avail- 
able. 

An analysis of the progenitor has to take into account 
the B-type companion discovered by Munari et al. H2[)()2ajl . 
It accounts for about half the brightness of the progeni- 
tor. It seems most reasonable to assume that V838 Mon 
and its B-type companion form a binary system. The 
main argument obviously comes from the observed po- 
sitions. From the instrumental crosses of stars seen on 
the HST images taken in September-December 2002 
I http: / /hubblesite.org/ newscenter/archive/2003 /10/| see 
also Bond et al. "20 035 O'^^ can deduce that the central 
object in the B images (dominated by the B-type com- 
panion) very well coincides with the central object in the 
I image (dominated by V838 Mon itself) and that both 
stars cannot be separated by more than ~ O'.'l. In the HST 
field (83" x 83") there are ~10 field stars of similar bright- 
ness as V838 Mon before outburst and the B-type com- 
panion. In this case the probability that due to a random 
coincidence one of these stars is separated by ^ O'.'l from 
V838 Mon is ^ 10~^. Next, as discussed below, the bi- 
nary hypothesis leads to a consistent interpretation of the 
observational data of the progenitor. Observational deter- 
mination of the distance and reddening to V838 Mon itself 
and its B-type companion, summarized and discussed in 
Tylenda (2005), give consistent results, in the sense that 
there is no significant difference in the results for both ob- 
jects. Therefore in most of our discussion we assume that 
V838 Mon and its B-type companion are at the same dis- 
tance and suffer from the same interstellar extinction. In 
some cases, however, we relax this assumption and discuss 
the consequences of that. 

As mentioned above, from spectroscopy Munari et al. 
Il2002afl have identified the hot companion of V838 Mon 
as a typical B3 V star. Indeed their photometric results 
obtained in September and October 2002, i.e. V = 16.05, 
B -V ^ 0.68 and [/ - B = -0.06 (Munari et al. 
can be well reconciled with the standard B3 V colours 
(see Schmidt-Kaler ITI)82|I provided that the object is red- 
dened with Eb-v = 0.90 (see open symbols and crosses 
in Fig. nil). 

The brightness of the B-type companion can also be 
deduced from the photometric results of Crause et al. 
(P005). Between September 2002 and January 2003 (AJD 
529-668 in Table 2 of Crause et al.) V838 Mon was prac- 
tically constant in U BV while steadily brightening in 
R and /. This behaviour of the object was also noted 
by Munari et al. (|2002al |2005') and indicates, in accord 
with their spectroscopic results, that the UBV magni- 
tudes were dominated by the B-type companion during 
this time period. From the most reliable results of Crause 
et al., i.e. those for dates not marked with an asterisk 
in their Table 2, obtained between AJD 529-668 one de- 
rives (mean value ± standard deviation) V = 16.26±0.02, 
B-V ^ 0.51 ± 0.02, and [/ - B ~ -0.1. Thus the object 
is by ^0.2 mag. fainter in V than in Munari et al., while 
the B — V value if interpreted with the B3 V standard 
gives i?_B_y = 0.72. In this case however the object seems 



In the present paper we collect and discuss the data 
available on the progenitor of V838 Mon. This includes 
the archival photometric measurements done in the op- 
tical and infrared before 2002, results of analysis of the 
evolution of the light echo after the eruption, as well as 
available data on regions of interstellar matter (ISM) near 
the position of V838 Mon. In the case of erupting stars, 
conclusions drawn from the progenitor usually are very im- 
portant for constraining the mechanism of the eruption. 
An analysis of the observational data for V838 Mon during 
and after its eruption is done in another paper (Tylenda 

2. The analysis of the photometric data 

Table^lists the photometric results for V838 Mon prior to 
its outburst. Columns (1) and (2) give the names and the 
effective wavelengths of the photometric bands. The mag- 
nitudes and the error estimates in column (3) are from 
the references given in the last column. Optical magni- 
tudes have been taken from Kimeswenger et al. ( 2002 ) and 
Goranskij et al. ( 2004) . Munari et al. (2002b. .2005) have 
also estimated magnitudes of the V838 Mon progenitor. 
However, the results given in these two references differ 
by 1 mag. We do not take them into account as it is 
not clear what caused such large differences (Munari et al. 
I2()05l do not comment on this). However, if one takes mean 
values from Munari et al. I)2002bl l2QQ^ they do not sig- 
nificantly differ from those quoted in Tabled The object 
has also been observed in infrared surveys. JHK magni- 
tudes can be found in the 2MASS data while the DENIS 
experiment measured the UK bands. 

Note that different measurements have been based 
on observations taken at different epochs. However, the 
fairly constant B magnitude obtained in Goranskij et al. 
H200 4') between 1928-1994 shows that the progenitor of 
V838 Mon was not significantly variable. 

As can be seen from Table Q for four photometric 
bands we have two independent measurements. In the case 
of the B and J magnitudes the agreement is good. The 
values in the / and K bands are discrepant by ~0.2 mag- 
nitude. As it is difficult to judge which result is more reli- 
able, for futher analysis we have adopted mean values in 



R. Tylenda et al.: On the progenitor of V838 Monocerotis 



3 



13 



14 



15 



16 



17 
-0.5 



1 

- B3V + 


B1.5V 




1 


- Eg-V = 


0.90 




















X 

< 


- 




X 






1 




(a) - 



0.5 



log A(^tlTl) 



13 



14 



n3 



15 



16 



17 
-0.5 



1 

- B4V + 


A0.5V 




1 


- Eg-V = 


0.71 






- 






- 


- 




X 


X 

X _ 




X 

X 








1 




(b) - 



0.5 



log A(/j,m) 



Fig. 1. Spectrophotometry of the progenitor of V838 Mon. Full symbols - observed magnitudes from Tabled Part (a) 
(left panel): UBV magnitudes of the B-type companion (open symbols) and Eb-v = 0.90 adopted from Munari et al. 
H2Ul)2a|l . Crosses - standard B3 V star, full curve - standard B3 V and B1.5 V stars co-added. Part (b) (right panel): 
UBV magnitudes of the B-type companion (open symbols) and Eb-v — 0.71 derived from Crause et al. (|2Uy5| (see 
text). Crosses - standard B4 V star, full curve - standard B4 V and AO. 5 V stars co-added. General note: all the 
model spectra (crosses and full curves) have been reddened with the respective values of Eb-v- 



to be not blue enough in C/ — S. A better agreement with 
the above results derived from Crause et al. is obtained 
for a B4 V standard and Eb-v = 0.71 (see open symbols 
and crosses in Fig.^). 

From the beginning of October 2002 V838 Mon has 
also been measured by Goranskij et al. f2004) . From their 
results obtained in October-December 2002 one derives 
V = 16.17 ± 0.05 and B - y = 0.64 ± 0.13 (no U mea- 
surements have been done during this time period) . Thus 
the V magnitude is in between the values of Munari et al. 
and Crause et al., while the B — V value is closer to that 
of Munari et al. and implies Eb-v = 0.84 ± 0.14. A large 
scatter in the B measurements of Goranskij et al. should 
however be noted. 

Later on in this section we consider two cases depend- 
ing on whether the photometric data for the B-type com- 
panion are adopted from Munari et al. (2005) or from 
Crause et al. I)2()()5|l . The differences in the magnitudes 
between these two references are extreme (the data from 
Goranskij et al. 120041 are in between them) so these two 
cases allow us to see how the results of our analysis de- 
pend on uncertainties in the photometry of the B-type 
companion. 

Figure Q presents our interpretation of the available 
photometric data done assuming that both V838 Mon and 
the B-type companion have the same reddening. In the 
discussion we also assume that both components are at the 
same distance, namely that they form a binary system. In 
both parts of the figure full symbols display the observed 



magnitudes from Tabled The best fits, shown with full 
curves, have been made using the least square method and 
the intrinsic photometric colours for the main sequence 
stars taken from Schmidt-Kaler l|1982|l . Johnson ( 195^. 
Koornneef I^^M^ and Bessell & Brett l(Tn55|l (for more 
details on the fitting procedure see Tvlenda 12005(1 . 

In part (a) open symbols show the UBV photometry 
of the B-type companion taken from Munari et al. H2()()5|l 
fitted with a standard B3 V star shown with crosses. The 
full curve presents the best fit to the full points obtained 
with a standard B1.5 V star added to the B3 V compan- 
ion. Both spectral components have been reddened with 
Eb^v = 0.9. The ratio of the luminosity of the B1.5 com- 
ponent to that of B3 is 1.9. This ratio is somewhat too 
low for the B1.5 V and B3 V stars but given the uncer- 
tainties in the observational data we can conclude from 
Fig. n(a) as follows. The progenitor of V838 Mon was a 
binary system consisting of two early B main sequence 
stars. V838 Mon itself was probably somewhat brighter, 
hotter and more massive than its companion. The system 
is young, i.e. ^ 2 x 10^ yrs (main sequence lifetime of a 
9 Mq star, typical for B2 V). Note that it is excluded 
that V838 Mon was an evolved B1.5 star as then it would 
have been significantly more luminous than the B3 main 
sequence companion. 

Fig.^b) adopts the parameters of the B-type compan- 
ion derived from the photometry of Crause et al. H2005|l . 
i.e. a B4 V star reddened with Eb-v = 0.71. In this case 
in order to reproduce the observational data for the pro- 
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genitor an AO. 5 V standard star has been added to the 
B4 V companion. The luminosity ratio of the A compo- 
nent to the B one is 0.43. This is much larger than the 
ratio for the main sequence of the same types which, ac- 
cording to Schmidt-Kaler H1982|l . is ~ 0.02.^ Also the pos- 
sibility that the A component was an evolved star, e.g. a 
giant evolving towards the red giant branch (RGB), the 
AGB, or a post-AGB star, can be ruled out. In this case 
it would be expected to have been initially (while being 
on the main sequence) more massive and thus, at present, 
significantly more luminous than the B4 main sequence 
companion. Therefore the only possibility within the bi- 
nary hypothesis is that the AO. 5 star is in the pre-main- 
sequence phase. The system would thus be very young. 
Judging from the luminosity of the A- type component, 
~ 550 Lq if 1300 Lq is assumed for the B4 V compo- 
nent, its mass would be ~ 5 Mq and the age of the system 
would be of - 3 X 10^ yrs (Ihen HM^ . 

In summary, although the uncertainties in the observa- 
tional data for the progenitor and for the B-type compan- 
ion do not allow us to unambiguously identify the nature 
of V838 Mon the above discussion allows us to put rather 
narrow constrains if the most probable hypothesis of bina- 
rity is adopted. In this case V838 Mon is a system consist- 
ing of two intermediate mass stars. V838 Mon itself cer- 
tainly was not an evolved star, e.g. RGB, AGB, post-AGB. 
It is either slightly more massive than its B-type compan- 
ion, i.e. 8-10 Mq, and was on the main sequence prior 
to eruption, or is somewhat less massive, ^ 5 Mq, be- 
ing in the pre-main-sequence phase. The system is young, 
with the age estimated between 3 x 10^ and 2 x 10^ yrs. 
The abundances in V838 Mon obtained by Kipper et al. 
H2004|l are reminiscent of those in the so-called HAEBE 
stars (e.g. Acke & Waelkens 12004(1 which are believed to 
be more massive analogues of the T Tauri stars. Therefore 
it is likely that the V838 Mon system is still partly em- 
bedded in the interstellar complex from which it has been 
formed. Indeed, as discussed in Sect. 13.31 near the position 
of V838 Mon there are several star-forming regions with 
radial velocities close to that of V838 Mon and its B-type 
companion. This also fits well the conclusion of Tylenda 
H2004|l and Sect. l3.ll that the circumstellar dust producing 
the light echo of V838 Mon is most probably of interstellar 
origin. 

Munari et al. ( 2002|) have made an analysis of the pho- 
tometric data for the V838 Mon progenitor similar to ours. 
Their conclusion is qualitatively similar to ours in the 
sense that the progenitor was an early-type star. However, 
contrary to our main-sequence or pre-main-sequence hy- 
pothesis, Munari et al. conclude that the V838 Mon 
outburst was that of an evolved star of initial mass of 
^ 65 Mq, at present in a region occupied by Wolf-Rayet 

^ From this result one may consider that V838 Mon was an 
A-type main sequence star prior to eruption, thus rejecting 
the binarity hypothesis. However in the case of the AO. 5 main 
sequence star it would be at a distance of 3.0 kpc which is too 
low given the distance estimates from the light echo. 



stars in the HR diagram and having ~ 50 000 K. 
However, in a case like this we should see a bright HII 
region surrounding V838 Mon. The observed light echo 
(discussed in Sect. I3.l|l shows that there is a lot of diffuse 
matter extending from 0.1 pc up to at least ~ 4 pc. 
The 50 000 K star of Munari et al. (2005), assuming 
Eb-v = 0-9 and a distance of 8 kpc, would have a lumi- 
nosity of 3 X 10** Lq. Using model results of Stasihska 
(1990), for abundances depleted by a factor of 2 relative 
to standard values (V838 Mon lies at the outskirts of the 
Galactic disc) , we can estimate that a star like this would 
be able to ionize the surrounding matter up to Rg — 8 pc 
if its density is nn = 10 H atoms cm~^ (i?s scales as 
njj^^"^). The emission line spectrum would be dominated 
by [OIII] and Balmer lines ([OIII]A5007A/H/3 ~ 8), while 
the H/? luminosity would be ^ 230 Lq. For an ob- 
server (at 8 kpc and Eb-v = 0.9) it would look like 
a nebula with a diameter of 7' and an H/3 flux of 
~ 6xl0~^^ ergs~^ cm~^. The resultant H/3 surface bright- 
ness of ^ 4.5 X 10^^^ erg s~^ cm~^ arcsec^^ is typical for 
many extended planetary nebulae, e.g. those in the Abell 
II1966|I catalogue (observed H/3 fluxes and nebular diame- 
ters can be found in Acker et al. 11992(1 . Thus the nebula 
would be rather easy to discern observationally, especially 
that in Ha and [OIIIJA5007A it would 8 times brighter 
than in H/3. Yet no emission-line nebula have been discov- 
ered around the position of V838 Mon (Orio et al. 120021 
Munari et al. l2002b(l . Thus the idea of Munari et al. pUUSIl 
that V838 Mon prior outburst could have been as hot as 
50 000 K is not consistent with the observations. From the 
observed lack of any significant emission nebula around 
V838 Mon we can conclude that before the outburst the 
star was cooler than ^ 30 000 K, i.e. of a spectral type not 
earlier than BO. Munari et al. 1(2005(1 also consider that the 
progenitor could have had T^ff 25 000 K (although they 
argue that this is not likely). This solution is practically 
the same as our case of a B1.5 star in Fig.^a) which we 
interpret as an early B-type main sequence star. 

As discussed above it is evident that V838 Mon was 
not a typical red giant nor an AGB star prior to eruption 
if V838 Mon and its B-type companion form a binary sys- 
tem. The only way to reconcile the RGB or AGB hypothe- 
sis is to assume that the B-type companion has nothing to 
do with V838 Mon and that the coincidence of the two ob- 
jects in the sky is purely accidental. Then one may assume 
that V838 Mon is less reddened than the B-type compan- 
ion and a cooler star can be fitted to the observations. Let 
us consider that the B-type companion has the parame- 
ters derived from the observations of Crause et al. 1(2005(1 , 
i.e. B4 V reddened with Eb~v = 0.71, as then the fits 
give later spectral types for V838 Mon than if the results 
of Munari et al. l(20()2a(l were adopted. Let us also use 
the standard supergiant spectra (intrinsic colours taken 
from the same references as the main sequence ones) to 
model the contribution from V838 Mon. This is more rel- 
evant with the RGB/ AGB hypothesis and also results in 
later spectral types from the fits than the main sequence 
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spectra. Eb-v — 0.5 seems to be a lower limit for the 
extinction towards V838 Mon (see discussion of different 
observational determination in Tvlenda l2U(J5|l . Assuming 
this value, the best fit to the observations is obtained for 
the spectral type Fl (effective temperature ~ 7500 K). If, 
in spite of observational determination, the extinction is 
pushed to its hmit, i.e. Eb-v = 0.0 is assumed, the fit 
gives G7 (effective temperature ^ 4700 K). Thus there 
is no way to reconcile an M-type star with the observa- 
tional data. This conclusion is obvious if one realises that 
the B — R colour for an unreddened M-type star is > 2.7 
while that of the V838 Mon progenitor was ~ 1.3. 

From the above results we can firmly conclude that 
V838 Mon was not an AGB star. If one still does not want 
to leave the AGB hypothesis and argues that it best ex- 
plains the existence of the circumstellar matter seen in the 
echo and infrared images, then the only way to reconcile 
it with the photometric data is to say that V838 Mon had 
quite recently left the AGB and prior to eruption was in 
the post-AGB phase. However in this case its luminosity 
would be close to 5 x IO^Xq and its distance would have to 
be ~ 55 kpc in the case of the Fl type and Eb-v = 0.5 or 
even ~ 90 kpc if one prefers G7 and Eb-v = 0.0. Thus a 
typical spectral types of the post-AGB stars (F-G) would 
require rather unacceptable conditions, i.e. a very low red- 
dening and a very large distance putting V838 Mon in the 
extreme outskirts of the Galaxy. For the observationally 
acceptable range of the extinction, i.e. Eb-v between 0.7 
and 1.0 (Tylenda 12003, V838 Mon would have been of 
the B9-B0 type and, assuming the typical post-AGB lu- 
minosity as above, its distance would be between 7 and 
30 kpc. Thus the only possibility within the AGB - post- 
AGB hypothesis, which is not excluded by the photomet- 
ric data, is that V838 Mon was a B-type post-AGB star. 
However, as discussed above, the binarity with the B-type 
main sequence component is, in this case, excluded in spite 
of similar estimated distance ranges (7.5-12.5 kpc for the 
B-type companion, see Tvlenda l2UU5fl . close values of in- 
terstellar extinction and the same positions in the sky of 
both objects. If one adds that the duration of the phase 
when the post-AGB star can be classified as B-type is typ- 
ically lO'^ years it is clear that this solution is extremely 
improbable. 

Finally let us discuss the giant hypothesis. As dis- 
cussed above it is certain that V838 Mon was not a typi- 
cal RGB star, i.e. of K-M spectral type, prior to eruption. 
The latest acceptable spectral type, obtained assuming the 
lower limit of Eb-v = 0.5, is F0-A5. At the lower limit of 
the distance of 5 kpc f Tvlenda [2D05;1 the star would have 
a luminosity of 40-50 Lq which is more or less consis- 
tent with the standard giant luminosities for these spec- 
tral types fSchmidt-Kaler il982l . For the more probable 
reddening, i.e. Eb-v — 0.7 — 1.0, we have to move to the 
B types and correspondingly larger luminosities and dis- 
tances. Thus the hypothesis that V838 Mon was a giant 
before its eruption is not excluded provided that it was an 
early (A-B) type giant. The binarity with the B-type com- 
panion is excluded in this case (as discussed above). The 



star would be quite massive, ^ 2.5Mq, not far evolved 
from the main sequence and thus be in a fast evolutionary 
phase (time scale ^ 5 x 10^ years). The object would be 
quite rare in the stellar population although not as rare 
as the B-type post-AGB case considered above. The dis- 
cussed case would however have difficulties in explaining 
the origin of the circumstellar matter seen in the light 
echo. The wind from an early type giant would not be 
enough. On the other hand, the giant, being significantly 
older (most probably at least as old as 10^ years) than the 
B-type binary system considered above, would have little 
chance to still reside in a dense interstellar cloud. 



3. Circumstellar and interstellar enviroment 

3.1. The light echo 

The phenomenon of light echo observed in V838 Mon dur- 
ing and after the outburst suggests that there is much 
dusty matter in the vicinity of the object. Since the light 
echo works as a sort of scanner, an analysis of the light 
echo images in different epochs should provide detailed in- 
formation on the dust distribution near the object which 
would be important for constraining the nature of the ob- 
ject. Unfortunately, in spite of numerous images obtained 
at different observatories (including HST) no elaborate 
study of the light echo has been done as yet. So far the 
most detailed, but still very simple, analysis has been done 
by Tylenda H2004|l on five echo images obtained with HST 
between 30 April and 17 December 2002. His conclusion 
is that the dust distribution does not show any signs of 
spherical symmetry and that dust is likely to be of in- 
terstellar origin rather than due to past mass loss from 
V838 Mon. 

We have extended the analysis of Tylenda H2004|l 
using two recent echo images obtained on 21 Oct. 2003 at 
the USNO (http://www.ghg.net/akelly/v8381ar3.jpa) 

a nd on 8 Feb. 2004 with the HST 

( http : / /hubblesit e . org / newscenter / newsdesk / archive / releases /2004/ 10, 

On these images we have measured the outer rim of the 

light echo. Then a least square fit of a circle to the 

measurements has been done in the same way as in 

Tylenda ('2001. The results of the fits are given in the 

last two lines of Table |21 The first five lines in the table 

repeat the results from Table 1 of Tylenda f2004) as the 

uncertainties there have been slightly overestimated. The 

first column of Tabic [21 shows the time of observations, 

^obs, given in days since 1 January 2002. The radius of 

the echo, 6, and its uncertainty are given in the second 

column. The next two columns show the {x, y) position of 

the centre of the fitted circle relative to the central star. 

Note that x points to west while y is to north. The last 

column gives the angular distance of the echo centre from 

the central star. All the results are in arcsec. Following 

Tylenda (2004) we adopt in our analysis that the zero 

age of the echo is — 34 days (since 1 Jan. 2002). 

Filled symbols in the left panel of Fig. [21 display the 
evolution of the echo radius, 9, with time. The full curve 
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time (days since 1.01.2002) distance (kpc) 

Fig. 2. Left panel: The best fit of a plane model to the observed evolution of the light echo radius. Symbols with error 
bars - the values and uncertainties of 9 from column (2) of Table|5] Full curve - model predictions for d = 11.1 kpc and 
zq = 7.1 pc (see text). Right pannel: The 95% confidence region (hatched) of the plane model fitted to the observed 
light echo expansion (symbols in the left panel). 



Table 2. Results of fitting a circle to the outer edge of 
the light echo of V838 Mon. Time of observations, iobs, is 
in days since 1 January 2002. Results are in arcsec. 



^obs 



yc 



120.0 


18.55±0.62 


-0.83±0.90 


0.46±0.86 


0.95±0.88 


140.0 


20.89±0.61 


-1.12±0.88 


0.50±0.84 


1.23±0.86 


245.0 


30.32±0.83 


-2.27±1.21 


0.37±1.15 


2.30±1.18 


301.0 


33.68±0.84 


-2.36±1.23 


1.11±1.17 


2.60±1.20 


351.0 


36.66±0.90 


-2.61±1.31 


1.58±1.23 


3.05±1.27 


659.0 


52.10±1.15 


-6.05±1.65 


3.12±1.61 


6.80±1.63 


769.0 


56.67±1.19 


-5.91±1.66 


5.37±1.71 


7.99±1.69 



shows the best fit to the data of a plane slab model, i.e. 
of Eq. (17) in Tylenda piHljl . obtained for d ~ 11.1 kpc 
and Zq ~ 7.1 pc, where d is the distance between the light 
source and the observer while zq is the distance of the 
dust slab from the source. However, similarly to Tylenda 
H2004|l . the minimum of the fit is quite shallow and 
extended along zq ~ • The right panel of Fig. 13 shows 
the 95% confidence region of the fit. From this figure one 
can conclude that the distance to V838 Mon is 5 kpc. 
Recently, Crause et al. (|2UU5|I have analysed the light 
echo evolution from their observations done at the SAAO. 
Their results obtained for the sheet model are well within 
the hatched region in the right panel of Fig. [3 

As can be seen from Table |21 the centre of the light 
echo has been migrating from the central object. This mi- 
gration, displayed in Fig.|21 has kept the same pattern for 
2 years. The following two conclusions can be drawn 
from Fig. 13 First, since the appearance of the light echo 



its center has been moving away from the central object 
in roughly the same (north-east) direction, as can be seen 
from the left panel of Fig. |3 Second, the distance of the 
echo centre from the central star has been increasing lin- 
early with time, as shown in the right panel of Fig. |3 

As discussed in Tylenda (, .2004.) , the fact that the light 
echo has an outer edge means that the dusty medium pro- 
ducing the echo has a boundary in front of the central star. 
However, the fact that the echo edge is not centered on the 
star and that the distance of the echo centre from the star 
does increase with time shows that this dust boundary 
is not spherically symmetric with respect to the central 
object. 

Following the theory of the light echo, as e.g. summa- 
rized in Tylenda (2004) , the only reasonable interpretation 
of the observed evolution of the outer edge of the light echo 
is that the dust boundary in front of V838 Mon is more 
or less in the form of a plane inclined to the line of sight. 
A linear fit to the observed evolution of the distance of 
the echo centre from the central object, i.e. last column in 
Table 13 or symbols in the right panel of Fig. gives the 
relation 

9c = (O'.'OlOe ± O'.'OOBO) (—^) (1) 

Iday 

where t ~ t^hs ^ to is the time since the zero age of the 
echo. This relation is shown by a full line in the right panel 
of Fig. 13 Assuming a distance of 8 kpc Eq. (Q, together 
with Eq. (7) of Tylenda H2004|l . implies that the normal 
to the dust surface is inclined to the line of sight at an 
angle of ^ 26°. This surface is at a distance of ~ 3.5 pc 
from V838 Mon (along the line of sight) and its portion 
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Fig. 3. Migration of the echo centre from the central star. Left panel: {x,y) positions of the echo centre from columns 
(3) and (4) of Table [5] The size of the symbols is proportional to the echo radius given in column (2) of Table |2| 
Note that the x and y axes point to west and north, respectively. The central star is at (x — 0,y = 0). Right panel: 
evolution of the distance of the echo centre from the central star, 0c, with time. Symbols - the data from column (5) 
in Table 13 Full line - a linear fit to the data as given in Eq. 



so far (i.e. till Feb. 2004) illuminated by the light echo 
has dimensions of ~ 4.4 x 4.9 pc. It is difficult to imagine 
that such a large flat surface of the dusty medium could 
have been produced by mass loss from V838 Mon itself. 
Instead, it suggests that the dusty medium in the vicinity 
of V838 Mon is much more extended than the illuminated 
part, thus most probably being of interstellar origin. 

In April 2002, when V838 Mon faded in the optical, the 
light echo started developing an asymmetric hole in the 
centre. As analyzed in Tylenda 1)2004(1 this clearly shows 
that there is a dust-free region around V838 Mon and that 
this empty region is strongly asymmetric. The inner edge 
of the dusty region in the southern directions is at 0.10- 
0.15 pc from the central object whereas in the opposite 
directions it is at least 10 times further away. It can be 
noted that the IRAS PSC source 07015-0346 coincides 
with the position of V838 Mon within a position elipse of 
48" X 10" (see also Kimeswenger et al. I2()02|l . The mea- 
sured fluxes at 100 fim and 60 /im are 4.6 Jy and 1.4 Jy, 
respectively, while at 25 /im and 12 fim the catalogue gives 
only an upper limit of 0.25 Jy. When fitted with a sim- 
ple dust emission model, i.e. emissivity proportional to 
X~^Bx{Td), the IRAS fluxes give a dust temperature, Tj, 
of ^ 30 K. For a central source of 10^ Lq (two B3 V stars, 
see Sect. ^ this value of dust temperature is reached at 
~ 0.2 pc (see e.g. Eq. 7.56 in Olofsson EUOl- Thus the 
IRAS fluxes can be consistently interpreted as due to in- 
ner parts of the dusty region inferred from the light echo 
analysis. In particular, they give evidence that there is 



no significant amount of dust at distances ^ 0.1 pc, thus 
confirming the existence of the central dust-free region. 

The strongly asymmetric central dust-free region 
would be very difficult to understand if the echoing dust 
were produced by a past mass loss from V838 Mon. The 
hole would imply that mass loss stopped a certain time 
ago, e.g. 10"* yrs for the 0.10-0.15 pc inner dust rim if a 
wind velocity of 10 km s^^ is assumed. However the hole 
asymmetry would imply that in the opposite direction ei- 
ther mass loss stopped 10 times earlier or the wind veloc- 
ity was 10 times higher. Neither of these two possibilities 
seems to be likely. 

Instead, as discussed in Tylenda H2004(l . the asymmet- 
ric hole is easy to understand if the echoing dust is of inter- 
stellar origin. Then it is natural to suppose that V838 Mon 
is moving against the ISM. If possessing a fast wind it 
would create a hole largely asymmetric along the direc- 
tion of the movement. Indeed, the structure of the inner 
edge of dust in Fig. 5 of Tylenda (2004) well resembles 
stellar wind bow shocks investigated in e.g. Van Buren & 
McCray l(TU55[| and Wilkin ijlH^ . The nearest rim in the 
southern directions, being at 0.10-0.15 pc from the central 
star, would correspond to a region where the stellar wind 
collides head-on with the ambient medium. 

Let us assume that a star losing mass at a rate, M^, 
and a velocity, v^, is moving in the ISM of number den- 
sity, hq, with a relative velocity of u,. A swept-up shell is 
created in the form of a bow shock and in the up-stream 
direction this takes place where the wind ram pressure 
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is comparable to that of the ambient medium (see e.g. 
Van Buren & McCrayCnHi Wilkin IBMIl . i.e. 



no ^ Ww '^w 



(2) 



where is the number density in the wind and is related 
to the mass loss rate in a standard way 



(3) 



where mn is the H atom mass while /u is the mean molec- 
ular weight in units of mn. Then Eq. |(2Jl yields a standoff 
distance, ro, from the star, i.e. 



ro ~ 0.15 



1/2 



10-9 Mq yr-1 
no 



VlOkms-i/ Vlcni-3/ 



1000 km s-i 

-1/2 



1/2 



pc. 



(4) 



.10 km s-i 

where = 1.4 has been assumed. 

A B3 main sequence star has a luminosity of 5.0 x 
10^ Lq (Schmidt-KalerCnHS). Thus, according to the re- 

the expected mass loss 



lation of Howard & Prinja (jE 
rate would be 7 x 10"^° Mq yr'^ If V838 Mon is a young 
binary system, as discussed in Sect.|2 its velocity relative 
to the ISM should be rather low. Assuming between 
3 — 10 km s^^ and the above mass loss rate, Eq. yields 
no ~ 1 — 10 cni^-^ if ro = 0.10 — 0.15 pc. This result is 
uncertain but it indicates that V838 Mon is imbedded in 
a relatively dense ISM. In principle, the density could be 
estimated from the echo brightness. Unfortunately there 
is no such estimate available, although the presence of 
the bright echo suggests that the density of the ambient 
medium must be significant. 

A fast stellar wind colliding with a circumstellar 
medium should produce X-rays. V838 Mon was observed 
with Chandra a year after the outburst by Orio et al. 
H2()()3|l . The object was not detected and the upper limit 
to the X-ray luminosity is ~ 0.13 Lq (for a distance of 
8 kpc). The kinetic power of a wind with parameters as 
above (7 x 10"^" Mq yr"^ expanding at 1000 km s^^) 
is ~ 0.06 Lq, thus it is below the X-ray limit, although 
only by a factor of 2. However, the X-ray luminosity from 
a colliding wind is usually much below the wind kinetic 
power (e.g. Soker & Kastner 2003 ). A hot bubble is usu- 
ally formed and, if the radiative cooling time is long, the 
shocked gas cools adiabatically by expansion of the bub- 
ble. In our case the cooling time is estimated to be above 
10* years, i.e. much longer that the life time of a B3 main 
sequence star. In the case of a bow-like inner edge of the 
circimstellar matter, as discussed above, a hot bubble need 
not be formed. Instead the shocked gas may flow along the 
edge and escape through the open side, thus expanding, 
cooling adiabatically and radiating very little X-rays. 

van Loon et al. (2004) have also analyzed the expan- 
sion of the light echo. Their analysis has been based on 
their measurements of the echo diameter on images from 
different sources. However, as also noted in Crause et al. 
I|2()()5|l . their diameters are systematically smaller by fac- 
tor of 2.5-2.7 than any other measurements available in 



the hterature (Munari et al. l2002b[ Tvlenda OUMI Crause 
et al. 20(131 see also numerous individual measurements 
in lAU Circ. in 2002), including our present results in 
Table 12 It is curious that these authors do not note this 
discrepancy and do not comment on it. In any case it can 
be concluded that the whole analysis of the light echo 
made by van Loon et al. is questionable as it has been 
based on wrong data. In particular, their rather spec- 
ulative interpretation of the outer edge of the echo in 
October 2003 - February 2004 as being produced by scat- 
tering under right angles does not hold. With the correct 
values of the echo diameter in these dates this interpreta- 
tion would imply a distance of ~ 2 kpc (and not 5.5 kpc 
as written in van Loon et al) . which is much too low com- 
pared with any other distance estimates from the light 
echo evolution (Bond et al. Tylenda I^HMl Crause 

et al. I 



3.2. Infrared and CO shells of van Loon et al. 



van Loon et al. (2004, hereafter lvLER,S|l from their analy- 
sis of the IRAS and MSX images and the CO maps claim 
that V838 Mon is surrounded by three shells. The inner- 
most, seen from the MSX, is highly irregular and has di- 
mensions of ~ The second one would be the elliptical 
one referred from the IRAS with dimensions of 15 — 20'. 
The largest one, having a diameter of ~ 1°, is suggested 
from the CO maps. On this basis fyTERSI conclude that 
V838 Mon is a low mass AGB star experiencing thermal 
pulses. As we have shown in Sect. El the photometric data 
on the V838 Mon progenitor exclude the AGB hypothesis. 
However, even if the latter is not taken into account, we 
find severe problems with the results of IvLERHI and their 
interpretation. 

First questions arise when analyzing the innermost 
structure seen in the band A (8.3 /im) of MSX. Unlike 
the two outer shells (IRAS and CO) showing an ellipti- 
cal or circular symmetry this structure is very distorted 
and does not show any kind of symmetry with respect to 
V838 Mon. It does not resemble objects involving AGB, 
like planetary nebulae, envelopes of AGB and post- AGB 
stars. Usually, as in planetary nebulae for example, the 
ISM affects external regions so one would expect to see 
more distortion in the IRAS and CO shells than in the 
MSX structure. 

The dimensions of the IRAS and CO shells seem to 
be too large to be compatible with the hypothesis that 
they have been produced by an AGB mass loss. Adopting 
a distance to V838 Mon of 8 kpc fTvlenda l^UUljl the ra- 
dius of the IRAS shell is ~ 20 pc while that of the CO 
shell is ^ 70 pc. (Munari et al. 120051 argue for 10 kpc as 
the most probable distance, which would make radii even 
larger thus strengthening our conclusions below.) The 
largest observed AGB dust shells have radii of 2-3 pc 
(Speck et al. 12000(1 while the CO shells are usually well 
below 1 pc (Olofsson'20Q4j). We can thus conclude that the 
IRAS and CO shells, if real, are not typical AGB shells. 
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Below we show that it is unlikely that AGB shells could 
survive and be observed as fairly symmetric structures at 
distances of 20-70 pc from the central star. As a mass- 
losing star moves relative to the ISM, the shell's segment 
in the up-stream direction (the side facing the ISM) is 
slowed down, until it is stopped at a time ^stop, when the 
leading edge is at a distance rgtop from the star. Next the 
up-stream segment of the shell is pushed by the ISM to- 
ward the central star. At the same time the shell's segment 
in the down-stream direction is expanding at a constant, 
undisturbed rate. Soker et al. H1991fl have derived simple 
analytical expressions for these parameters which can be 
used here. 

Let a shell of mass, Mj, expand with a velocity Wcxp, 
and let be the relative velocity of the mass-losing star 
and the ISM. Also let pQ be the mass density of the ISM, 
and uq the total number density of the ISM. For a dis- 
tance of ~150 pc from the galactic plane, we can scale the 
ISM density with po = 10~^^gcm~^, which corresponds 
to riQ — O.lcm^^. Following Soker et al. H1991|l we define 
the radius of a sphere which contains an ISM mass equal 
to the shell mass 




The stopping distance of the up-stream shell's segment is 
given by (Eq. 5 of Soker et al. [TMT|l 

rstop =i?o [2a (l + a)]-i/3, (6) 

where a = v^,/vcxp- The time the up-stream segment 
reaches this maximum distance is 



= J^a-'/\ (7) 

^cxp 

As the down-stream segment expands undisturbed the di- 
ameter of the shell along the stream at tstop is 

4top = Ro [(2a (1 + a))-i/3 + ^-2/3] 

After reaching its maximum distance from the central star 
at Tstop , the up-stream segment of the shell is pushed by 
the ISM toward the central star. 

For an AGB star Vcxp — 10 km s^^. The typical 
star-ISM velocity at ~ 150 pc from the galactic plane is 

> 10 km s~^. For a = 1 and assuming Ms — 0.1 Mq, 
we get from Eqs. (0, ® and ||SJ) rgtop — 1.6 pc and 
c'stop = 4 pc. Even for an extreme case of no = 0.01 cm^^, 
Mg = 1 Mq (which is a generous upper limit for an AGB 
shell, especially if it were ejected in a thermal pulse from a 
low mass star, as suggested in lvLER^ and a = 0.25 (say, 
expansion velocity of 20 km s,~^ and — 5 km s~^), we 
find Tstop = 13.5 pc and dstop = 43 pc. Given the observed 
diameter of the CO shell of ~140 pc we can conclude that 
an AGB shell would have been seriously disturbed by the 
ISM before reaching these dimensions, its up-stream part 



would have to be significantly brighter (because of sig- 
nificant accretion of the matter from the ISM) than the 
opposite part and the central star would very likely be 
now observed outside the up-stream rim. All this is not 
observed. 

The above estimates are supported by observations of 
planetary nebulae. From Eq. jSJ we see that when a typ- 
ical planetary nebula shell of 0.2 Mq reaches a radius of 
~ 3 pc, it is expected to be highly distorted by the ISM. 
Indeed, examining the list of planetary nebulae interacting 
with the ISM compiled by Tweedy & Kwitter ifl^Ml their 
Table 3) , we find that all the planetary nebulae in this list 
have radii < 5 pc. Most of the large planetary nebulae are 
highly distorted; not is only the central star not at the cen- 
ter, but the shells are neither circular nor elliptical. Some 
planetary nebulae, like NGC6826 NGC 2899 and A 58 
(surrounding the final helium shell flash star V605 Aql), 
have very large, diameters ~ 10 — 40 pc, IRAS structures 
around them (Weinberger & Aryal 120041 Clayton & De 
Marco [T997|l . They are all severely distorted. Clayton & 
De Marco H1997|l argue that structures of this size are 
swept-up ISM dust, rather than AGB mass-loss shells. 

The above analysis rises a question: are the shells 
claimed in IvLER'^ indeed real and related to V838 Mon? 

It is difficult to discuss the nature of the emission seen 
in the MSX image as it has been recorded only in the A 
band (8.3 /im) image. No emission is seen in the bands C 
(12.1 /im), D (14.6 urn) and E (21.4 ^m). This is perhaps 
due to the highest sensitivity of band A. 

The elliptical structure around V838 Mon in the IRAS 
image shown in Fig. 1 of vLERS. at first sight looks con- 
vincing. However, in the whole field of this figure it is easy 
to fit several ellipses of similar sizes and similar orienta- 
tions as the one drawn bv IvLERSi One possible explana- 
tion is that the image pattern seen in Fig. 1 of IvLEEItI 
might be spurious, i.e. of instrumental and/or image pro- 
cessing origin. Another likely interpretation is that this is 
a general pattern of the interstellar diffuse emission in this 
region and thus it has nothing to do with V838 Mon. The 
discussed region is a part of an extended infrared emission 
related to several molecular clouds and HII regions near 
the direction to V838 Mon (see Sect. l3.3|l . Whether or not 
a part of this emission is physically related to V838 Mon is 
an important question but it cannot be decided just from 
the image. 

Fig. 3 of lvLEH'^Tl has been derived from a compilation 
of surveys in the CO 1-0 line made by Dame et al. H2001|l . 
However, if one takes the composite map from Fig. 2 of 
Dame et al. and expands near the position of V838 Mon 
the resultant image is not the same as that in lvLERSI In 
particular, there is no emission to the left and upper-left 
of the position of V838 Mon [i.e. I > Z(V838 Mon) and 
b ^ 6(V838 Mon)], so no bubble-like structure is seen. 
Apparently Dame et al. f20011 considered this emission as 
statistically insignificant. Indeed, the data for this region 
come from a low resolution (0?5) and low sensitivity survey 
of Dame et al. (TMf). Thus the emissions seen in Fig. 3 of 
IvLERSI and coming from this last survey is uncertain and 
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Fig. 4. The 100 fim IRAS image centered on the position 
of V838 Mon (north is up, east to the left). The grid shows 
the galactic (/, b) coordinates. Numbers show the positions 
of the ISM regions from Table El 



might be spurious. This is supported by the fact that two 
large, faint patches seen to the left and upper-left of the 
V838 Mon bubble in Fig. 3 of vLERS (not seen in Fig. 2 
of Dame et al. I2001|) cannot be identified with any known 
CO cloud while two known CO regions, namely regions 
(9) and (13) in TableEKsee Sect. l3.3|l . are not seen in the 
image of lvT.FRSI 

3.3. Interstellar medium 

The presence of the light echo proves that there is dusty 
matter around V838 Mon. As argued in Tylenda ( 2004) 
and in Sect. IXTl of the present paper this matter is likely 
to be of interstellar character. This notion is supported 
by the conclusion of Sect. |21that V838 Mon is likely to 
be a young binary system, as well as by the fact that 
having the Galactic coordinates, I = 2 17? 80, 5 ~ +1?05, 
the object is located near the Galactic plane. Therefore it 
is important to investigate observational data on the ISM 
in the vicinity of V838 Mon. In this section wc discuss the 
available data from the IRAS and CO surveys. 



Figure 0| shows the IRAS image at 100 fim centered 
at the position of V838 Mon. At this wavelength prac- 
tically all emission seen comes from dust (interstellar or 
circumstellar). As can be seen from the figure V838 Mon 
is located in a faint diffuse emission probably related to 
bright regions near the Galactic plane. The numbers in 
the figure show the positions of molecular and HII regions 
listed in Table |31 which are within ^ 1?5 of the position 
of V838 Mon. The table gives the galactic coordinates of 
the regions, their usual names, values of taken from 
Wouterloot & Brand (5,989j, Vlsr resulting from CO line 
observations, and types of the regions. The references to 
the data are given in the last column of the table. 

As can be seen from Table 01 the regions marked in 
Fig. 21 can be divided into two groups from the point of 
view of their positions, brightness and I^sr. The brightest 
regions near and slightly below the Galactic plane (regions 
5, 6, 8, 12) have Vlsr in the range of 20-30 km s'^. The 
heliocentric radial velocity of V838 Mon is not well known 
as it has been estimated from outburst spectra. The results 
are within 55-65 km (Kolev et al. 120021 Kipper at al. 
I2004L M. Mikolajewski - private communication). That of 
the B-type companion is ~ 64 km s^^ (T. Tomov - private 
comminication) . Using the results of Dehnen & Binney 
l|1998(l this can be transformed to Vlsr = 44-54 km s"^ 
Thus the above regions are most probably at much smaller 
distances (2-3 kpc if interpreted with the Galactic rotation 
curve of Brand & Blitz [T^^ than V838 Mon and they are 
simply seen in front of the object. 

However the regions lying above the Galactic plane (re- 
gions 1, 2, 3, 4, 7, 9, 13) have Vlsr between 47-57 km s~^. 
When interpreted with the rotation curve of Brand & Blitz 
(jl993|l their distances are in the range of 6-8 kpc. Thus 
these regions are located much closer to V838 Mon and 
a physical relation between one of them and the matter 
seen in the light echo is quite possible. 

Region (10), whose apparent position is closest to that 
of V838 Mon, has been listed in Magakian L2003j) as a re- 
flection nebula related to a 9 magnitude B9 star HD 53135 
(LS —03 15) estimated to be at a distance of ~ 2 kpc 
(VogtCnZl Kaltcheva & Hilditch lMIUIl . Thus this region 
is probably located well in front of V838 Mon. 

Interstellar Na I lines in the spectrum of V838 Mon 
during eruption showed two components at heliocentric 
velocities of ~ 37 and ^ 64 km (Zwitter & Munari 
IMI^ Kolev et al. IMM Kipper et al. lM)^ . When trans- 
formed to Vlsr the figures become ^ 26 and ^ 53 km s~^. 
Thus both line components can be interpreted as due 
to ISM related to the above two groups of the ISM re- 
gions. Detailed observations of the vicinity of V838 Mon 
in molecular lines might be important for discussing the 
nature of V838 Mon. 



4. Discussion and summary 

The goal of this paper is to use available observational 
data on the progenitor and enviroment of V838 Mon to 
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Table 3. Interstellar regions within ~ 1?5 from the position of V838 Mon 



No. I b Name Vlsr Notes References 



(°) (°) (K) (km/s) 



1 


216.5+0.5 






47. 


MC 


MAE 


2 


216.5 +1.2 






48. 


MC 


MAB 


3 


217.0+0.9 


FT 84 




51. 


MC 


Av, DHT 


4 


217.2+0.5 




2.6 


52. 


MC 


Av, WB, MAB 


5 


217.3+0.0 






27. 


MC 


MAB 


6 


217.4-0.1 


FT 87, BFS57 


18.9 


26. 


MC 


Av, WB 


7 


217.4+0.3 


FT 88, BFS58 


3.7 


50. 


MC 


Av, WB, MAB 


8 


217.6-0.2 


FT 89, BFS59 


10.4 


26. 


MC 


Av, WB 


9 


217.6+2.4 




6.3 


55. 


MC 


Av, WB 


10 


217.9+0.9 








RN 


Ma 


11 


218.0+0.2 






23. 


MC 


MAB 


12 


218.1-0.3 


S287, FT 91 


18.2 


26.-30. 


HIT, MC 


Av, WB, MAB 


13 


218.7+1.8 


IC466, S288 


6.3 


57. 


HIT, MC 


Av, WB 



Notes: HII - HII region, MC - molecular cloud, RN - reflection nebula. 

References: Av - Avedisova (I2UU2II . DHT - Vlsr estimated from original data of Dame et al. II2UU1II . Ma - Magakian II2UU3II . 
MAB - May et al. llTWft . WB - Wouterloot & Brand llBH^ . 



better constraint the nature of its eruption. Below we sum- 
marize and discuss our main findings and conclusions. 

(1) The nature of the progenitor. In Sec. El we have 
analyzed the photometric data available for the progen- 
itor. Most likely the progenitor was a young binary sys- 
tem consisting of two intermediate mass (5-10 Mq) stars. 
V838 Mon itself was either a main sequence star of similar 
mass as its B-type companion or a slightly less massive 
pre-main-sequence star. The system is very wide as the 
B-type companion observed today does not seem to be af- 
fected by the eruption. From the maximum photospheric 
radius of V838 Mon during eruption (Tylenda 2005| we 
can estimate that the separation of the components is 
> 3 X 10'^ Rq so the orbital period is > 12 years. The 
B-type companion was probably not involved in the erup- 
tion of V838 Mon, at least directly. The hypothesis of a 
young binary system is also supported by the position of 
the object near the Galactic plane and the conclusion of 
Sect.Elthat V838 Mon is probably embedded in the ISM. 

A less likely hypothesis is that the presently observed 
B-type companion does not form a binary system with 
V838 Mon. In this case, other than a B-type main sequence 
star, the progenitor could have been an A-B spectral type 
giant evolving from the main sequence or a B-type post- 
AGB star. These two possibilities however involve a short 
(giant in the Hertzsprung gap) or very short (post-AGB) 
evolutionary phase. As it is an old object in these cases it 
would not be expected to reside inside or close to dense 
ISM regions. 

We can safely excluded the possibility that before erup- 
tion V838 Mon was of spectral type K-M, so it could not 
have been a typical RGB or AGB star. 

(2) The light echo and the Galactic enviroment of 
V838 Mon. In several studies the light echo was used 
to argue that the light-reflecting dust was expelled by 
V838 Mon in previous eruptions (e.g Bond et al. i2003i 



lvLERS|l . As argued by Tylenda J20 045, and discussed here 
in Sect. 13. II the data strongly suggests that dust is of ISM 
origin. The dust structure derived from the echo analy- 
sis does not show any hint of spherical symmetry. On the 
contrary, the outer boundary of the echoing dust in front 
of the object can be approximated by a plane at a dis- 
tance of ^ 3.5 pc from V838 Mon and inclined at an 
angle of ^ 26° to the line of sight. The strongly asym- 
metric dust-free region in the near vicinity of V838 Mon, 
inferred from the central hole in the echo, is interpreted as 
produced by the V838 Mon progenitor (and possibly its 
B-type companion) moving relative to the local ISM and 
sweeping out the medium by its fast wind. As discussed in 
Sect. 13.31 there are several interstellar molecular regions 
seen in the IRAS image and CO surveys probably located 
near V838 Mon. The local ISM seen in the light echo is 
therefore likely to be related to one or some of them. 
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